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We evaluate the dressed Polyakov loop for hot quark matter in strong magnetic field. To compute 
the finite temperature efi^ective potential, we use the Polyakov extended Nambu-Jona Lasinio model 
with eight-quark interactions taken into account. The bare quark mass is adjusted in order to 
reproduce the physical value of the vacuum pion mass. Our results show that the dressed Polyakov 
loop is very sensitive to the strenght of the magnetic field, and it is capable to capture both the 
deconfinement crossover and the chiral crossover. Besides, we compute self-consistently the phase 
diagram of the model. We find a tiny split of the two aforementioned crossovers as the strength of 
the magnetic field is increased. Concretely, for the largest value of magnetic field investigated here, 
eB = 19m^, the split is of the order of 10%. A qualitative comparison with other effective models 
and recent Lattice results is also performed. 
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I. INTRODUCTION 

The nature of the Quantum Chroniodynamics (QCD) 
vacuum is one of the most intriguing aspects of modern 
physics. Besides, it is very hard to get a full understand- 
ing of its properties, because its most important char- 
acteristics, namely chiral symmetry breaking and color 
confinement, have a non-pcrturbative origin, and the use 
of perturbative methods is useless. One of the best strate- 
gies to overcome this problem is offered by Lattice QCD 
simulations at zero chemical potential (see [l|-[j] for sev- 
eral examples and see also references therein). At van- 
ishing quark chemical potential, it is almost established 
that two crossovers take place at nearly the same tem- 
perature; one for quark deconfinement, and another one 
for the (approximate) restoration of chiral symmetry. It 
is still under debate whether two crossovers should oc- 
cur at exactly the same temperature, see for example the 
report in Ref. [2|. 

An alternative approach to the physics of strong in- 
teractions, which is capable to capture some of the 
non-pcrturbative properties of the QCD vacuum, at the 
same time being easy to manage mathematically, is the 
Nambu-Jona Lasinio (NJL) model Q, see also Refs. [g 
for reviews. In this model, the QCD gluon-mediated in- 
teractions are replaced by effective interactions among 
quarks, which are built in order to respect the global 
symmetries of QCD. Since dynamical gluons are absent 
in this model, it is not a gauge theory. However, it shares 
the global symmetries of the QCD action; moreover, the 
parameters of the NJL model are fixed to reproduce some 
phenomenological quantity of the QCD vacuum: in its 
simplest version, the pion decay constant, the vacuum 



raoul. gatto@unige.chl 



,ruggieri@yukawa.kyoto-u.ac.jp| 



pion mass and the vacuum chiral condensate are repro- 
duced. Therefore, it is reasonable that the main charac- 
teristics of its phase diagram represent, at least qualita- 
tively, those of QCD. 

Critically speaking, the worst aspect of the NJL model 
is that it lacks confinement: massive quark poles of the 
quark propagator are present at any temperature and/or 
chemical potential. It is well known that color con- 
finement can be described in terms of the center sym- 
metry of the color gauge group and of the Polyakov 
loop [3| , which is an order parameter for the center sym- 
metry. Motivated by this property, the Polyakov ex- 
tended Nambu-Jona Lasinio model (P-NJL model) has 
been introduced [1, Q , in which the concept of statisti- 
cal confinement replaces that of the true confinement of 
QCD, and an effective interaction among the chiral con- 
densate and the Polyakov loop is achieved by a covariant 
coupling of quarks with a background temporal gluon 
field. In the literature, there are several studies about 
various aspects of the P-NJL model. Its phase structure 
with two fiavors and symmetric quark matter has been 
investigated in Refs. |10l - [l3| : PNJL model with a Van der 
Monde term has been considered in [1J| ; phase structure 
with 2-1-1 fiavors has been studied in Refs. [1^; possi- 
ble realization of the quarkyonic phase [la] has been dis- 
cussed in |15l.ll7l|: mass dependence of the phase diagram, 
and a possible emergence of the quarkyonic phase, is in- 
vestigated in [l8|; phase diagram with imaginary chemi- 
cal potential has been studied i n [l9l . l2g | ; dual quark con- 
densate has been computed in [2l|; neutral phases have 
been investigated in |22| ; phase diagram with asymmetric 
quark matter have been studied in [23| ; non-local exten- 
sion has been introduced in [24| ; role of eight-quark inter- 
actions in the PNJL context has been elucidated in j25| . 

The modification of the QCD vacuum, and of its ther- 
mal excitations as well, under the influence of external 
fields, is an attractive topic. Firstly, it is extremely inter- 
esting to understand how an external field can modify the 



main characteristics of confinement and spontaneous chi- 
ral symmetry breaking. Lattice studies on the response 
to external magnetic fields can be found in |26l - [28l |. QCD 
in chromo-magnetic fields has been investigated on the 
Lattice in [2^[3y|. Previous studies of QCD in magnetic 
fields, and of QCD-like theories as well, can be found 
in Refs. [3ll - J34l |. A self-consistent model calculations of 
magnetic catalysis and of deconfinement pseudo-critical 
ternperature in magnetic field, has been performed firstly 
in [35| within the PNJL model, and then in [3a| using the 
Polyakov extended quark-meson model. Effective models 
in chromo-magnetic fields have been considered in [37[. 
Besides, strong magnetic fields might be produced in non- 
central heavy ion collisions [3a l39| . In this case, it has 
been argued that the non-trivial topological structure of 
thermal QCD gives rise to Chiral Magnetic Effect (CME) 

Beside the Polyakov loop, it has been suggested [43 
that another observable which is an order parameter 
for the center symmetry, hence for confinement, is the 
dressed Polyakov loop. From the mathematical point of 
view, the dressed Polyakov loop is built from the canoni- 
cal (called thin) Polyakov loop, by dressing it with higher 
order loops, which wind once around the compact tem- 
poral direction. In this context, the order of a loop is 
given by its length; the thin Polyakov loop corresponds to 
the shortest one. The dressing becomes more important 
when quark masses are finite (Polyakov loop is an ex- 
act order parameter for confinement-deconfinement only 
in the ideal case of static quarks with infinite masses). 
In Refs. [ij] the dressed Polyakov loop has been com- 
puted within the scheme of truncated Schwinger-Dyson 
equations, with a model for resummed quark-gluon ver- 
tex and in-medium gluon propagator computed on the 
Lattice. Within the Nambu- Jona Lasinio model, in which 
the QCD interaction among quarks is replaced by a con- 
tact four-fermion interaction. Si has been computed at 
finite temperature and chemical potential in [4J] . Finally, 
the dressed Polyakov loop has been computed within the 
PNJL model in [2l| at finite temperature. 

In this article, we compute the phase structure and the 
dressed Polyakov loop of hot quark matter at zero chem- 
ical potential, in an external magnetic field. To compute 
the effective potential, we rely on the PNJL model of 
strongly interacting quarks. We will limit ourselves to the 
one- loop approximation (saddle point), which is enough 
to draw a phase structure. The novelty of the present 
article is many fold. Firstly, we introduce the eight-quark 
interaction |45l-l48l| in the PNJL model in an external 
magnetic field (previous studies of the PNJL model in 
magnetic and chromomagnetic fields neglected this kind 
of interaction). Within the Nambu- Jona Lasinio model, 
it has been shown that the eight-quark interactions natu- 
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rally lowers the pseudo-critical temperature for (approxi- 
mate) chiral symmetry restoration. Magnetic catalysis in 
the NJL model with multi-quark interactions has been in- 
vestigated in [i^. However, in those studies, the compu- 
tation of quantities relevant for deconfinement crossover 
is lacking. On the other hand, in [23 the Polyakov ex- 
tended NJL model with multi-quark interaction has been 
investigated, but without magnetic field. It is of interest, 
then, to study the response of quark matter to magnetic 
fields in the framework of the PNJL model with eight- 
quark interaction. In doing this, we will consider quarks 
with finite values of bare mass, fixed to reproduce the 
vacuum pion mass, while in a previous study [S^ this 
problem was studied only in the chiral limit. 

Moreover, we compute the dressed Polyakov loop, Ei, 
in a magnetic field. Along this line, we anticipate one 
of our results, namely that the dressed Polyakov loop, 
Si, is capable to feel both the Polyakov loop and the 
chiral condensate crossovers, whatever the strength of 
the magnetic field is. This occurs despite the tiny split 
of the two crossovers, which we observe at sufficiently 
strong magnetic field strength. Therefore, in view of an 
effective theory for finite temperature QCD in terms of 
just one order parameter, our results are encouraging. 

The plan of the paper is as follows. In Section II, we 
present the model we use. In Section III, we show and 
discuss our numerical results. Finally, in Section IV, we 
draw our conclusions. 



II. DRESSED POLYAKOV LOOP IN THE 
EFFECTIVE MODEL 

In this article, we model quark matter by the following 
Lagrangian density 

C = qii'y'^Df, - mo) q + g^ [(qqf + {qhsTqf] 

+98 [{qqf + iqil5Tqf]'' , (1) 

which corresponds to the NJL lagrangian with multi- 
quark interactions |45| . The covariant derivative embeds 
the quark coupling to the external magnetic field and 
to the background gluon field as well, as we will see ex- 
plicitly below. In Eq. ([1]), q represents a quark field in 
the fundamental representation of color and flavor (in- 
dices are suppressed for notational simplicity); mo is the 
bare quark mass, which is fixed to reproduce the pion 
mass in the vacuum, ttItt — 139 McV. Our interaction in 
Eq. ([T]) consists of a four-quark term, whose coupling g^ 
has inverse mass dimension two, and an eight-quark term, 
whose coupling constant gs has inverse mass dimension 
eight. 



The evaluation of the bulk thermodynamic quantities requires we compute the quantum effective action of the model. 
This cannot be done exactly. Hence, we rely ourselves to the one-loop approximation for the partition function, which 
amounts to take the classical contribution plus the fermion determinant. The one- loop thermodynamic potential of 



quark matter in external fields has been discussed in I35|,|37[, in the case of canonical antiperiodic boundary conditions; 
following |2l| . it is easy to generalize it to the more general case of twisted boundary conditions: 
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In the previous equation, a = ga-iqq) = '2ga{uu); fc is a 
non-negative integer which labels the Landau level; ak = 
Sko+2{l — Sko) counts the degeneracy of the /c— th Landau 
level. We have put 

UkiPzf^pl + 2\qfeB\k + M^, (3) 

with M ~ nio — 2(7 — Aa^gs/g^ ■ The arguments of the 
thermal exponentials are defined as 



£± =ujk{pz) ± 



i{ip — tt) 



(4) 



with if defined in Eq. ^. 

The vacuum part of the thermodynamic potential, 
Cl(T = 0), is ultraviolet divergent. This divergence is 
transmitted to the self-consistent equations which deter- 
mine the chiral condensate and the expectation value of 
the Polyakov loop. In this article, we use a smooth reg- 
ularization procedure by introducing a form factor (?a(p) 
in the diverging zero-point energy. Our choice of ^a ip) is 
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(5) 



we choose two values of N, namely N = 5 and N ~ 7. 

The potential term U[P,P,T] in Eq. © is buih by 
hand in order to reproduce the pure gluonic lattice 
data [lO[. Among several different potential choices [43 
we adopt the following logarithmic form [3, [l3 ; 



U[P,P,T] =T^ 



a{T) 



$$ 



b{T) In 1 - QPP + A{P^ + P^) - 3(PP) 



(6) 



with three model parameters (one of four is constrained 
by the Stefan-Boltzmann limit). 



a{T) = ao + ai ( ^ j + as f Y 



h{T) = h,['^ 



(7) 



The standard choice of the parameters reads [1^ ; 

ao = 3.51, ai = -2.47, as = 15.2 , 63 = -1.75. 

(8) 
The parameter Tq in Eq. ([6]) sets the deconfinement scale 
in the pure gauge theory, i.e. Tc = 270 MeV. 

Following [421 j i'^ order to define the dressed Polyakov 
loop, we work in a finite Euclidean volume with tem- 
perature extension /3 = 1/T. We take twisted fermion 
boundary conditions along the compact temporal direc- 
tion. 



q{x,p) = e-''Pq{x,Q) , ^€[0,2^] 



(9) 



while for spatial directions the usual periodic boundary 
condition is taken. The canonical antiperiodic boundary 
condition for the quantization of fermions at finite tem- 
perature, is obtained by taking </? = tt in the previous 
equation. The dual quark condensate, E„, is defined as 



S„(to, F) = 



27r 



dip e^"^ 



■{qq)c 



(10) 



where n is an integer. The expectation value {•)g denotes 
the path integral over gauge field configurations. An im- 
portant point is that in the computation of the expec- 
tation value, the twisted boundary conditions acts only 
on the fermion determinant; the gauge fields are taken 
to be quantized with the canonical periodic boundary 
condition. 

Using a lattice regularization, it has been shown in [42| 
that Eq. (jlOp can be expanded in terms of loops which 
wind n times along the compact time direction. In par- 
ticular, the case n = 1 is called the dressed Polyakov 
loop; it corresponds to a sum of loops winding just once 
along the time direction. These correspond to the thin 
Polyakov loop (the loop with shortest length) plus higher 
order loops, the order being proportional to the length 
of the loop. Each higher order loop is weighed by an in- 
verse power of the quark mass. Because of the weight, in 
the infinite quark mass limit only the thin Polyakov loop 
survives; for this reason, the dressed Polyakov loop can 
be viewed as a mathematical dressing of the thin loop, 
by virtue of longer loops, the latter being more and more 
important as the quark mass tends to smaller values. 



TABLE I. Parameters of the model for the two choices of the 
UV-regulator. 





A (MeV) 


mo (MeV) 


g^ (MeV) -2 


58 (MeV)-** 


N = 5 


588.657 


5.61 


5 X 10^^ 


6 X 10-^2 


N = 7 


603.475 


5.61 


4.92 X 10"'' 


6.8 X 10^22 



If we denote by z an element of the center of the color 
gauge group, then it is easy to show that S„ —5- z"S„. 
It then follows that, under the center of the symmetry 
group Z3, the dressed Polyakov loop is an order param- 
eter for the center symmetry, with the same transfor- 
mation rule of the thin Polyakov loop. Since the center 
symmetry is spontaneously broken in the deconfinement 
phase and restored in the confinement phase (7| (in pres- 
ence of dynamical quarks, it is only approximately re- 
stored), the dressed Polyakov loop can be regarded as 
an order parameter for the confincment-dcconfincmcnt 
transition as well. 

For later convenience, we scale the definition of the 
dressed Polyakov loop in Eq. ([TU| , and introduce 
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where (7[ip) corresponds to the expectation value of the a 
field computed keeping twisted boundary conditions for 
fermions. 



III. NUMERICAL RESULTS 

In this Section, we show our results. The main goal to 
achieve numerically is the solution of the gap equations. 



an 

9^ 



= 0, 



on 

dP 



= 



(12) 



This is done by using a globally convergent algorithm 
with backtrack [50] ■ From the very definition of the 
dressed Polyakov loop, Eq. ((TO)) . the twisted boundary 
condition, Eq. (jH]), must be imposed only in D^. There- 
fore, we firstly compute the expectation value of the 
Polyakov loop and to the chiral condensate, taking Lp = it. 
Then, in order to compute the dressed Polyakov loop us- 
ing Ec^. (|lip , we compute the <y9-dcpendcnt chiral conden- 
sate using the first of Eq. ([T^ . keeping the expectation 
value of the Polyakov loop fixed at its value eA, ip - 
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In this study, we report results obtained using the UV- 
rcgulator specified in Eq. ([5|) with N ~ 5 and N — 7. As 
expected, there is no qualitative difference among the 
pictures that the two regularization schemes lead to. As 
a consequence, concrete results are shown only for the 
case N = 5; for what concerns the case N = 7, we collect 
the pseudo-critical temperatures in Table [TTl We have 



also checked that the results are qualitatively unchanged 
if we use a hard cutoff scheme instead of the smooth UV- 
rcgulator. The parameter set for both cases is specified in 
Table m In the case N = 5, they are obtained by the re- 
quirements that the vacuum pion mass is mjr = 139 MeV, 
the pion decay constant /^r = 92.4 MeV and the vacuum 
chiral condensate (uu) ~ (—241 MeV)'^. In this case, the 
chiral and deconfinement pseudo-critical temperatures at 



zero magnetic field arc Tq 



rpP 



175 McV. Simi- 



larly, for the case N = 7, the chiral and deconfinement 
pseudo-critical temperatures at zero magnetic field are 
T^ = 176 MeV and T^ = 175 MeV, respectively; the 
zero temperature chiral condensate at zero magnetic field 
strength is fixed to (uu) « (-246 MeV)^. 

We remark that the main effect of the eight-quark in- 
teraction in Eq. ^ is to lower the pseudo-critical temper- 
ature of the crossovers. This has been already discussed 
several times in the literature [4^, |43 , in the context of 
both the NJL and the PNJL models. Therefore, it is not 
necessary to discuss it further here, while at the same 
time we prefer to stress the results that have not been 
discussed yet. 

In order to identify the pseudo-critical temperatures, 
we have define the effective susceptibilities as 



XA = (m^)s 



dA 



A = a,P, El 



(13) 



Strictly speaking, the quantities defined in the previous 
equation are not true susceptibilities. Nevertheless, they 
allow to represent faithfully the pseudo-critical region, 
that is, the range in temperature in which the various 
crossovers take place. Therefore, for our purposes it is 
enough to compute these quantities. In Equation (|13p . 
the appropriate power of ttt-tt is introduced just for a mat- 
ter of convenience, in order to have a dimensionless quan- 
tity; therefore, g = if A = cr. Si, and g = 1 ii A = P. 



A. Condensates and dressed Polyakov loop 

From now on, we fix A^ = 5 unless specified. The 
results for this case are collected in the form of three- 
dimensional plots in Fig. [T] (for the case N = 7 the plots 
do not differ qualitatively). In the left panel, we plot 
the chiral condensate (uu)^/'^, the expectation value of 
the Polyakov loop, and the dressed Polyakov loop Si, 
as a function of temperature and magnetic field. In the 
right panel, we show the contour plots of the raw data 
of the effective susceptibilities. The lighter the color, the 
higher the susceptibility. In the contour plots, the verti- 
cal axes correspond to temperature (measured in MeV); 
the horizontal axes represent the magnetic field eB/m^. 

We slice the three dimensional plots in Fig. [T] at fixed 
value of the magnetic field strength, and show the re- 
sults in Fig. [5J where we plot the chiral condensate 
S = |(uu)|^/^ (upper panel), the Polyakov loop (middle 
panel) and Si (lower panel) as a function of temperature, 
for several values of the applied magnetic field strength. 
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FIG. 1. Left panel. Chiral condensate, Polyakov loop and dressed Polyakov loop as a function of temperature and magnetic 
field, for the case A'' = 5. Right panel. Contour plots of the raw data of the effective susceptibilities. The lighter the color, the 
higher the susceptibility. Vertical axes correspond to temperature (in MeV); horizontal axes represent magnetic field eB /m\. 
For the dressed Polyakov loop susceptibility, the bifurcation of the critical region is evident. 
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FIG. 2. Left panel. Chiral condensate S — \{uu)\^'''^ (upper panel), Polyakov loop (middle panel) and Ei (lower panel) as 
a function of temperature, for several values of the applied magnetic field strength, measured in units of m^. In the figures, 
N = 5 corresponds to the order of the UV- regulator in Eq. (0. Right panel. Effective susceptibilities, defined in Eq. (|13p . as a 
function of temperature, for several values of eB. Conventions for lines are the same as in the left panel. 



measured in units of m^ . In the right panel, we plot fits 
of the effective susceptibilities in the critical regions, as 
a function of temperature. The fits are obtained from 
the raw data, using Breit-Wigner-likc fitting functions. 
The details of the fitting procedure are not relevant for 
the present discussion. For graphical reasons, in Fig. [T] 
we plot the chiral condensate with its sign; on the other 
hand, in Fig. [5] we take the absolute value of this quantity. 

The qualitative behavior of the chiral condensate, and 
of the Polyakov loop as well, is similar to that found in 
a previous study within the PNJL model in the chiral 



limit [35|. Quantitatively, the main difference with the 
case of the chiral limit, is that in the latter the chiral 
restoration at large temperature is a true second order 
phase transition (in other model calculations it has been 
reported that the phase transition might become of the 
first order at very large magnetic field strengths [34|). 
On the other hand, in the case under investigation, chiral 
symmetry is always broken explicitly because of the bare 
quark masses; as a consequence, the second order phase 
transition is replaced by a smooth crossover. 

Another interesting aspect, observed also in the chi- 



ral limit [35|, is that the Polyakov loop crossover tem- 
perature, is less sensitive to the strength of the mag- 
netic field than the same quantity computed for the chi- 
ral condensate. It is useful, for illustration purpose, to 
quantify the net shift of the pseudo-critical temperatures, 
for the largest value of magnetic field we have studied, 
eB = 19rn^. In this case, if we take TV = 5 (for N = 7 
the results are similar), then the two crossover occur si- 
multaneously at eB = 0, at the temperature Tq = T^ = 
175 MeV; for eB = IQm^, we find T^ = 219 MeV and 
Tp = 190 MeV. Therefore, the chiral crossover is shifted 
approximately by 25.1%, to be compared with the more 
modest shift of the Polyakov loop crossover, which is 
^ 8.6%. 

The split of the two critical temperatures at a so large 
value of the magnetic field strength is only of 15%; on the 
Lattice, no split is observed [2g|, and a modest increase of 
the critical temperature is measured. Therefore, we are 
in partial agreement with the Lattice results, in the sense 
that the raising of the critical lines is observed also in our 
model calculation; for what concerns the split of the two 
crossovers, we can take our 0(10%) split as a consequence 
of the crudeness of the model at hand. On the Lattice, 
the smaller pion mass used is of the order of 200 MeV [23 . 
We have verified that our qualitative picture is unchanged 
if we increase artificially the vacuum pion mass up to 
this value. In passing, we notice that using a running 
coupling as in |52| . but adding at the same time two 
further free parameters in the model, we expect a better 
agreement with the Lattice. The reason is that in |52| . 
the coupling g^ is a function of the Polyakov loop, and it 
decreases as P is increased. As a consequence, near the 
Polyakov loop crossover temperature, the strength of the 
interaction is lowered, and a partial suppression of the 
chiral condensate is expected. Quantitatively, it is not 
clear a priori if the suppression is enough to rejoin the 
two crossovers; only a detailed numerical study can give 
the answer. We leave this important investigation to a 
future study. 

The tiny decoupling of the two crossovers found within 
the PNJL model, both in the chiral limit [35| and in the 
case of physical pion mass considered here, is observed 
also within the Polyakov quark-meson model |36|, when 
in the latter the zero point energy is considered (if the 
vacuum energy is subtracted, then the Polyakov loop and 
the chiral crossovers occur always simultaneously, but the 
pseudo-critical temperature is a decreasing function of 
eB, which seems in disagreement with the recent Lattice 
results [2g ; see also [51| for a recent discussion of the role 
of the vacuum energy within the quark- meson model). 
Since the Polyakov loop is coupled to quarks in the same 
manner both in the PNJL and in the PQM model, the 
tiny split of the two crossovers as eB is increased does 
not appear as an artifact of the PNJL model; instead, it 
seems to be a consequence of the link among the chiral 
condensate and the Polyakov loop, which is common in 
the two kinds of models. 

In the lower panels of Figures. [T] and [21 we plot the 



dressed Polyakov loop as a function of temperature, for 
several values of eB. Our definition, Eq. (fTTj) . differs 
from the canonical one [421 for an overall factor, which 
gives mass dimension one to our Ei. For small values of 
eB/rn^, the behavior of Si as temperature is increased, 
is qualitatively similar to that at eB = 0, which has been 
discussed within effective models in [2l|, |4J] . In partic- 
ular, the dressed Polyakov loop is very small for tem- 
peratures below the pseudo-critical temperature of the 
simultaneous crossover. Then, it experiences a crossover 
in correspondence of the simultaneous Polyakov loop and 
chiral condensate crossovers. It eventually saturates at 
very large temperature (for example, in [2l| the satura- 
tion occurs at a temperature of the order of 0.4 GeV, in 
agreement with the results of [4J|). However, we do not 
push up our numerical calculation to such high tempera- 
ture, because we expect that the effective model in that 
case is well beyond its range of validity. 

As we increase the value of eB, as noticed previously, 
we observe a tiny splitting of the chiral and the Polyakov 
loop crossovers. Correspondingly, the qualitative behav- 
ior of the dressed Polyakov loop changes dramatically: 
the range of temperature in which the Si crossover takes 
place is enlarged, if compared to the thin temperature 
interval in which the crossover takes place at the lowest 
value of eB (compare the solid and the dotted lines in 
Fig. [21 as well as the the lower panel of Fig. [1]). 

The effective susceptibility, dSi/dT, plotted in the 
lower right panel of Fig. [21 is qualitatively very interest- 
ing. We observe a double peak structure, which we inter- 
pret as the fact that the dressed Polyakov loop is capable 
to feel (and hence, describe) both the crossovers. If we 
were to interpret Ei as the order parameter for deconfine- 
ment, and the temperature with the largest susceptibility 
with the crossover pseudo-critical temperature, then we 
obtain almost simultaneous crossover even for very large 
magnetic field. If this were the case, then the Polyakov 
loop computed within the PNJL model, should be inter- 
preted only as an indicator of statistical confinement, and 
the deconfinement would be described by Ei. Of course, 
this picture would not contradict the well established pic- 
ture at zero magnetic field [BiHJ- Indeed, in the case of 
small eB, we find simultaneous crossover of chiral con- 
densate, Polyakov loop and dressed Polyakov loop. In 
the latter case, it would be just a matter of taste which 
quantity one uses to identify the deconfinement crossover. 
Even if it is tempting to give this kind of interpretation, 
which would lead to simultaneous crossover also at finite 
eB, it is very hard to accept it without more convincing 
microscopic arguments. Therefore, in the prosecution 
of this work, we prefer to associate the deconfinement 
crossover to that of the Polyakov loop. Nevertheless, the 
dressed Polyakov loop is a new quantity which is interest- 
ing to compute. In particular, the double peak structure 
in the Ei effective susceptibility, which is produced if the 
magnetic field is strong enough, offers the evidence that 
the dressed Polyakov loop is intimately related to both 
chiral condensate and (thin) Polyakov loop, and it is ca- 



TABLE II. 1 


Coefficients of the fit function defined in Eq. IQ^ . 




a 


a 


To (MeV) £ 


T^,N = 5 


2.4 X 10"^ 


1.85 


175 0.21 


Tp, N = 5 


2.1 X 10"^ 


1.41 


175 0.08 


T^,N = 7 


7.8 X 10-^ 


1.29 


176 0.19 


Tp,N ^7 


3.9 X 10-3 


1.08 


176 0.01 



pable to capture both the crossovers. The bifurcation of 
the dressed Polyakov loop susceptibiUty is impressive in 
the lower right panel of Fig. [TJ 



B. Phase diagram in the eB — T plane 

In Figure[21 we collect our results on the pseudo-critical 
temperatures for chiral and Polyakov loop crossovers, in 
the form of phase diagrams in the eB — T plane. The 
dashed line denotes the Polyakov loop crossover, and 
the dot-dashed line corresponds to the chiral crossover. 
The shaded area is the region, in the eB — T plane, in 
which quark matter is not statistically confined, but chi- 
ral symmetry is still broken by the chiral condensate. 
Temperature on the vertical axes are measured in units 
of the pseudo-critical temperature at zero field, which is 
To = 175 McV. Wc fit our data on the pseudo-critical 
temperatures by the law 



-^ c 
To 



1 



eBy 
To) 



(14) 



where A = a, P. Numerical values of the coefficients 
in Eq. (|14p for the various observables are collected in 
Table |TT1 As an estimator of the goodness of the various 
fits, we report in Table HIl the percentage error defined as 



e = 100 X ^ 



fA{Xi 



(15) 



where the sum runs over the data, {xi,yi) corresponds 
to a couple in the set of the data (cBjTa), and fA{xi) 
denotes the numerical value of the fit function evaluated 
at the data eB. 

The picture discussed in the previous Section is made 
clear by the phase diagrams in Fig. [31 We measure an 
increase of both deconfinement and chiral crossovers; the 
tiny split of the two critical temperatures is of the order 
of 10% for the largest value of the magnetic field strength 
considered here. 

It is instructive to compare our results with those ob- 
tained in a different model. The shape of the phase di- 
agram drawn in Fig. [3] is similar to that drawn by the 
Polyakov extended quark-meson model, see e.g. Fig. 13 
of Ref . [33 ■ In that reference, an interpretation of the 
split in terms of the interplay among vacuum and ther- 
mal contribution, is given. We totally agree with those 



arguments, which are reproduced within the PNJL model 
as well, as the final results on critical temperatures show. 
In the case of the quark-meson model, however, the pic- 
ture can change even qualitatively, depending on the fate 
of vacuum energy contribution. If they are not included, 
then a simultaneous first order transition is observed at 
every value of eB (only if eB is very small the transition 
is a smooth crossover), and the deconfinement tempera- 
ture as a function of the magnetic field strength decreases. 
This picture confirms the scenario anticipated in a pre- 
vious work [3J] . In the case of the PNJL model, we can- 
not reproduce the latter scenario, because of a technical 
reason: indeed, in our case the vacuum contribution can- 
not be subtracted (as a matter of fact, we do not have 
a further effective potential term at zero temperature, 
which leads to spontaneous breaking of chiral symme- 
try when vacuum quark contributions arc subtracted). 
Therefore, we can limit ourselves only to a comparison 
with the quark-meson model with vacuum contributions 
taken into account. 



IV. CONCLUSIONS 

We have computed, for the first time in the literature, 
the dressed Polyakov loop for hot quark matter in exter- 
nal magnetic field. To compute the finite temperature 
effective potential in magnetic field, we have used the 
Polyakov extended Nambu-Jona Lasinio model, with a 
logarithm effective action for the Polyakov loop. In the 
quark sector, we have used both a four-quark and an 
eight-quark interactions. Bare quark masses are fixed to 
reproduce the physical value of the vacuum pion mass. 
This model allows to treat self-consistently both chiral 
symmetry breaking and (effective, or statistical) confine- 
ment. We improve the previous work [35| in three ways: 
we set the vacuum pion mass to its physical value; we in- 
troduce eight-quark interaction; finally, we compute the 
dressed Polyakov loop. 

Our results on the dressed Polyakov loop, Ei, in mag- 
netic field show that this quantity is capable to describe 
both Polyakov loop and chiral crossovers. This is re- 
sumed in the double peak structure of the effective sus- 
ceptibility dSi/dT, see FigsH] and El Moreover, we 
find that Ei is capable to feel both the Polyakov loop 
crossover and the chiral condensate crossover, and sug- 
gests itself as the the possibly unique order parameter of 
effective QCD. 

The results on the pseudo-critical temperatures as a 
function of eB are resumed in the phase diagrams in 

Fig. [31 These results were anticipated in a previous 

work [35| in which only the chiral limit was considered, 
and the eight quark interaction was neglected. Our re- 
sults agree qualitatively with those of Ref. [3g , in which 
a quark-meson model coupled to the Polyakov loop is 
considered. 

As improvement of our results, it would be interest- 
ing to consider the effects of non- locality 12J]. In that 




eB/m/ 



FIG. 3. Phase diagram of the PNJL model in magnetic field. Dashed line denotes the Polyakov loop crossover; dot-dashed line 
corresponds to the chiral crossover. The shaded area is the region, in the eB — T plane, in which quark matter is not statistically 
confined, but chiral symmetry is still broken by the chiral condensate. Temperatures on the vertical axes are measured in units 
of the pseudo-critical temperature at zero field, which is To = f 75 MeV. The analytic form of the lines corresponding to Tp 
and T^ is specified by Eq (|14|l . The UV-regulator is that corresponding to A*' = 5. 



case, however, the computation of the fermion spectrum 
in the magnetic field is not trivial because of the non-local 
structure of the action. Another interesting possibility is 
the use of Montccarlo methods to compute the PNJL 
partition function in magnetic field, going beyond the 
saddle approximation. Encouraging results along this re- 
search line in the context of the PNJL model have been 
reported in [53 . Therefore, it might be interesting to 
extend the computation of [53| to the case of quarks in 
external magnetic field. Even more, we expect that the 
running coupling introduced by the Kyushu's group [52| 
would help (at least partly) to get closer crossovers in 
magnetic field. A numerical investigation of this subject 
is left to a future study. Finally, the extension of our cal- 
culation to finite quark chemical potential, and to quark 



matter in external chromo-magnetic fields, the latter be- 
ing motivated by Lattice results [23, [SO] , would deserve 
further attention. 
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